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Abstract: To get the total cross section of one interaction from its amplitude M., one needs 
to integrate \Ai\ 2 over phase spaces of all out-going particles. Starting from this paper, we will 
propose a new method to perform such integrations, which is inspired by the reduced phase space 
integration of one-loop unitarity cut developed in the last few years. The new method reduces one 
constrained three-dimension momentum space integration to an one- dimensional integration, plus 
one possible Feynman parameter integration. There is no need to specify a reference framework in 
our calculation, since every step is manifestly Lorentz invariant by the new method. The current 
paper is the first paper of a series for the new method. Here we have exclusively focused on massless 
particles in 4D. There is no need to carve out a complicated integration region in the phase space 
for this particular simple case because the integration region is always simply [0, 1]. 
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1. Introduction 



In the last few years there were great progresses in the evaluation of one-loop amplitudes for general 
field theory (see [jy and references within). One of such progresses is the unitarity cut method, which 
was initiated in |2], [3| and then pushed by Witten's "twistor program" M. One key achievement 
along this line is the reduced unitarity phase space integration using "holomorphic anomaly" 
More accurately, the reduced unitarity phase space integration is given by 1 

J d 4 L5 + (L 2 )5 + ((L - Kf)G(L) . (1.1) 

With the first delta-function 5 + (L 2 ) we can reduce the measure J d A L5 + (L 2 ) into the integration 
with variables t, A, A, where t is a affine variable and A, A are spinor variables. With the second 
delta- function 5 + ((L — K) 2 ), we can integrate t out and get 

(X\d\) \x\dx\ /(A, A) . (1.2) 

Before the work of ||, the evaluation of remaining two-dimensional integration is a very hard task 
and has blocked practical applications of unitarity cut method. From || is that people realized that 
the remaining integrations over A, A can be obtained by reading out residues of corresponding poles. 
In another word, there is no need for integrations and everything is just algebraic manipulation. 
This method is usually referred to as "Spinor Integration Method" or the spinor method for short. 

Clearly, the success of the spinor integration method is due to the presence of two delta- 
functions. Using this technique, we are able to perform any unitarity phase space integration at 
one-loop in pure 4 dimension [7], Originally, spinor is tightly related to null momentum in pure 
4D. However, for practical applications, it will be useful to generalize it to general D-dimension as 
well as massive particles. This goal has been achieved late in || [Kj . 

To get the cross section, we need to integrate the physical phase space 2 of all (or some) out-going 
particles, 

/ f 

where 5 4 (Pi n — P ou t) is the energy- momentum conservation condition and c is a function depending 
on the in-coming particles (for example, c _1 = 4EaEb\va — vb\ in 2 — > n progresses). For simplicity, 
we omit c in this paper. 



1 The 5 + ( ) means the integrated region with Iq > 0. 

2 To make things clear, we use unitarity phase space and physical phase space to distinguish these two cases. 
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Usually, integrations as such are difficult to perform, especially when there are many out-going 
particles. For LHC experiments, channels with four or even five out-going particles play important 



role in search of new physics (see, for example, ]TT[ and references in this paper). Aiming at this 
task, we will try to do the physical phase space integration using the spinor integration method. 

The advantage of the spinor method is that it reduces the constrained three-dimension mo- 
mentum space integration to an one- dimensional integration, plus one possible Feynman parameter 
integration. All remaining integrations are scalar type, i.e., the integrand is manifestly Lorentz 
invariant, so there is no need to specify a reference framework in our calculation. Furthermore the 
integration region in the phase space also becomes simpler. 

The current paper is the first of a serys of work we will take to complete our new method. In 
this paper, we will focus on massless particles in 4D. For this case, with familiar infrared/collinear 
divergences for massless out-going particles, there will be infinity after the phase space integration. 
To get sensible physical quantities, some regularization is needed. Due to this difficulty, the result 
in this paper is still far from the practical applications. It is the basis for all late work. We will deal 
with massive particles in the second paper and general D-dimension integration in the third part. 

For massless particles in pure 4D as focused on in this paper, some simplifications happen. For 
example, the integration region in phase space will always be [0, 1] in our method. This will be 
modified to be nontrivial functions of mass and energy when the particles are massive as to be 
presented in second paper. 

The main aim of this paper is to laid out the framework of our new method, so most examples 
in this paper are not for real cross sections. These examples is to demonstrate the salient character 
of spinor integration, such as frame independence and simple integration region. 

The outline of this paper is the following. In section 2, we transfer the physical phase space 
integrations into the form of spinor integrations. We start with the case when there are two and 
three outgoing particles, then generalize to cases with n > 3 recursively. The recursive feature (not 
recursion relation) is one of advantages of the spinor integration method. 

In the following three sections, our method has been demonstrated with simple examples with 
two, three and four outgoing particles, respectively. There will also a brief discussion on the 
IR/collinear divergence problem related to massless particles in section 3. 

In section 6, we summarize our results along with general discussions and remarks. 

There are two appendixes. In Appendix A the evaluation of pure four particle phase space 
directly using momentum components has been given as to compare with spinor integration method 
presented in main text. In Appendix B we have explained the unfamiliar spinor integration method 
from the point of view of integration in complex plane. From this point, the correctness and the 
power of this method becomes obvious. 
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Il 

Figure 1: Representation of the cut integral. Left and right legs are on-shell. Internal lines represent the 
cut propagators. 

1.1 The reduced phase space integration of unitarity cut 

As we have emphasized above, the key gradient of spinor method is the existence of two delta- 
functions. There is one situation where these two delta-functions arise naturally: it is the reduced 
phase space integration of one-loop unitarity cut. In this subsection, we review how to perform the 
phase space integration of one-loop unitarity cut by spinor method. 

Unitarity Cut: Consider the unitarity cut in the (i, i + 1, . . . ,j — 1, j)-channel, as shown in 
fig. [l| The cut integral is 



dfiA^ih, i.i • 1 j- 1../. Z 2 )A trcc (H 2 ), j + 1, j + 2, . . . , % - 2, % - 1, Hi)), (1.4) 



where dp, = dHidH25 ( - + \ll)5^ + \l2)5 4 (li+l2 — Pij) is the Lorentz invariant phase space measure of two 
light-like vectors (/i, I2) constrained by the energy- momentum conservation, and P^ = pi + . . . + pj. 
Notice that the integrand is similar to \Ai\ 2 in eq ( |1.3| ). So we can make full use of this point in 
our calculations. 

Using the 5 4 ( ) function to perform the I2 integration gives 



J 2 ) 



CW,..,;-u = / rf 4 W (+) (^ (+) (('i -Pi: 

xA^ih, i.i ■ 1 ./- 1../. l 2 )A trcc (—l 2 ,j + 1, j + 2, . . . , i - 2, i - 1, -h)(l.h) 



Note that the Lorentz invariant measure of a null vector I can be represented as a measure over 
R+ x CP 1 x CP 1 . The contour of integration is a certain diagonal CP 1 . Explicitly, one writes 
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La = t\ a \ a , and then 3 



d 4 W (+) (/ 2 )(«)~ / tdt (A dX) [dX A](«), 

JO </A=A 

where (•) represents a generic integrand up to an overall numerical factor. So we have 4 



t dt 



(X dX) 



X dX 



5^ P*-t(X\R,\X 



G(X, X,t). 



;i.6) 



;i-7) 



G(A, A, t) arises from the product of the two tree-level amplitudes in (|1.4|). In order to get G(A, A, t) 
in actual calculations, we have to write expressions of the form (•, I2) or [•, I2} in terms of I = l\ {l\ 
is always substituted by / when there is no possibility of confusion). A systematic way of doing this 
is as follows: 

A similar identity is valid for [•, l 2 ]. The factors [l 2 h] and (I2 h) all pair up in the end allowing for 
the use of the vector form of l^. This is because the product of amplitudes must be invariant under 
the scaling zXi 2 and z~ Y Xi 2 . For ( p..7|) , integrating over t yields 



(X dX) 



X dX 



X I Pij I A 



G A, A, 



P 2 



A I -P^ I A 



;i.9) 



The requirement of degree zero in A implies that G in ( |1.9|) can be written as a sum of terms in the 
form 5 









A 


a 




[x\Qi 


X 



(1.10) 



where R{, Qj are functions of external momenta. 

Canonical Splitting: In order to calculate the integral efficiently, we can reduce the integrand 
by separating the denominator factors with A as much as possible, at the cost of more terms. When 



3 Here we have used the QCD convention so measure is (A dX) [dX X] instead of the twistor convention (A dX) [A dX] 



4 Where we have written (l\ — Pij) 2 — Pfj ~ t (^X\Pij\X , thus we have used the QCD convention instead of the 
twistor convention. 

5 It is worth to emphasize that this general form is only true for the one-loop calculation. When wc apply our 
method to cross section evaluation, we will have other forms of input. The way to deal with them generally is given 



in Appendix B. 
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there is a product [a X][b A] in the denominator, multiply both numerator and denominator by [a b}. 
Applying Schouten's identity yields 



[i j] [k I] = [i k] [j I] + [i I] [k j] 



in the numerator with another factor [c A] (which must exist by homogeneity when the degree of 
A in denominator is equal to or more than three). Thus we get two terms with [a A] or [b A] in 
the numerator, canceling one of the denominator factors. The result, in terms of A, is a denomi- 



nator of the form Y[ r \A|Q r |A 



A A 



can be treated by writing 



in every term. Similarly ^A|Q r | A 

X\Q r \X = Qr A , where A Qr = -(Q r )aatf- 

Using this procedure repeatedly and noticing the integrand has a degree of —2 in A, we end up 
with two kinds of possible integrals: 



(A dX) 



A dX 



nliWiA 



1 



A | Pij | A 



B 



(A dX) 


A dX 




(a|^|A 


(A|Q r |A 



;i.i2) 



Xb can be changed to 



X, 



B 



da: 



(X dX) 



X dX 



X\xPij + (1 — x)Q r \X 



2 ' 



;i.i3) 



which reduces to X4 up to a Feynman parametrization integration. On the other hand, X4 is just 
the Yl i yi component of an auxiliary integration 6 

(A dX) ' 



T 



X dX 



X\R\X 



A|P y -|A 



2+k 



R = S ^ ViRi 



;i.i4) 



Now all integrations are reduced into the one given by T aux . It can in turn be written as a total 
derivative by using 



AIPIA 



/ 



1 



iv A] 



V 



n 



+ l)(X\P\rj] ^ |p| 



A 



n+l 



;i.i5) 



/ 



Thus, the evaluation of formula Q1-14Q is transformed into the reading of residues of 

k+l 

_^ / \ 1)1 \ 



X\R\X 



(k + l) (X\R\P tJ \X} / Mp 



\j\X 



k+l 



3 There are other ways to do it directly without using the auxiliary integration. See reference [Q. 
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where we need to sum up two possible poles from the factor (A|P|Pjj|A). The way to sum up these 
two poles has been discussed in detail, for example, in eq (70), (74) of |12| . In fact, the result can 
be written down directly as 

-[(-2P- J R+ v / A) fc+1 -(-2P- J R- VA) k+1 ] 

F= (1-lOJ 

(k + l) v / A(2P 2 ) fc+1 

with A = (2P ■ R) 2 — AP 2 R 2 . It is easy to see that numerator must be of the form Y2 a V / A 2a+1 (— 2P- 
R) k ~ 2a . The result is a rational function with the highest power of R to be k, as required. 

Before we end this part, let us emphasize that above procedure works for general one-loop cal- 
culation, but when we try to do the cross section calculation, we may meet new kinds of singularities 
and we need to generalize above procedure. The generalization has been discussed in Appendix B in 



some details. The basis idea is still to find formula like the one in ( \L.l3j ) and then take the residues. 



2. Spinor method for the physical phase space 

In this section we will establish the general framework for the physical phase space integration using 
spinor method. We will focus on massless particles in pure 4D. For massive particles and in general 
D-dimension, it will be discussed subsequent work. 

Before going into detail, let us make a simple observation: for two out-going particles, the 
physical phase space integration is the same as the unitarity phase space integration, which can be 
done directly by the spinor integration method. Difficulties arise when n > 3. We will first deal 
with the case n = 3 and then generalize to arbitrary out-going particles. 

Notice that we can rewrite 

d 3 pf 



, E( j <ll> f» ■(!>})■ (2.1) 
Then when there are only two out-going particles, we have 

J] j d 4 Pf 5 + (p 2 f )5\P m - Pl - p 2 ) = I d 4 Pl S + (pi)5 + (( P i - Pin?) ■ (2.2) 

The two delta-functions for J d 4 pi is exactly what one needs for the unitarity phase space integration, 
thus we can use the spinor integration method to evaluate the cross section efficiently. 

2.1 The phase space integration with three out-going particles 

When n = 3 we have 

3 



/ d 4 L l 5+(^))5 4 (^-^L l )/(L 1 ,L 2 ,L 3 ) 

i=l J 
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= Jd^S+iLl) J d 4 L 2 5 + (L\)5 + ((K - L ± - L 2 ) 2 )f(L 1 , L 2 , K-L x - L 2 ) 

The integration over L 2 with two delta-functions can be done by the spinor integration method. 
However, for f d A L\ we have only one delta- function. The same difficulty also arises when n > 3: 
for the integration over each momentum, we need another delta function in order to make use of 
the spinor integration method. Then the problem is how to rewrite the integration measure in such 
a form where one more delta-function shows up. The solution of this problem is the application of 
a Faddeev-Popov like method. We will show how it works in detail soon. 

Having the solution for the n = 3 case, the solution for arbitrary n > 4 can be obtained by 
doing the integration recursively one by one. Keeping (n — l)-momenta fixed, we can integrate the 
n-th momentum. After that, we keep (n — 2)-momenta fixed and integrate the (n — l)-th momentum 
and so on until n — 3. 

Now let us do the case n = 3 explicitly. The integration over L 2 is trivial. We will denote the 
result after this integration as f(Li). Then 

d 4 L 1 5 + (Ll)f(L 1 ) (2.3) 

which is a typical phase space integration and we rewrite it as 

dH 



I 



{2n / + (nm- (2.4) 

To get the standard form suitable for spinor integration method, we need to insert another 
5-function. To do so, let us consider the kinematics in more detail. The total energy-momentum 
tensor of out-going particles is given by K, which has following properties: K > for positive 
energy component and K 2 > 0. With these in mind, let us consider following expression 

I, = J dz5((zK - if) = J dz5(z 2 K 2 - z(2K ■£)+£*) = J dz^ \ 2z ^ - ■ i)y (2 ' 5) 



where 

(2K ■ I) ± y/(2K ■ I) 2 - 4KH 2 



(2.6) 



2K 2 

Let's determine the integration region of z for the physical system. In the center-of-mass frame, 
K = (E, 0, 0, 0), the out-going momentum £ = (Ei, P, 0, 0) and E\ = P 2 . Conservation laws ensure 
that the energy-momentum of all other out-going particles to be K ot h er = (E — Ei, — P, 0, 0). For 
physical particles we have K 2 ther > 0, i.e., (E — Ei) 2 > P 2 , or E > Ei + \P\. Then we get 

2E(E 1 ±P) E 1 ±P , N 
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When m = 0, we have E\ = \P\ and especially z_ = 0. Note that z = is always a solution when 
we combine the massless condition £ 2 = in (|2.5|) . Thus to have a physical meaningful solution, we 
can keep only z + solution. It is easy to see that 



< z , < 1 



(2.8) 



That is, I z defined in ( |2.5| ) makes sense physically only in region z G [0, 1]. With this consideration, 
I, can be defined 



dz 



8{z - z + ) 



1 



\y/(2K -£) 2 -4K 2 £ 2 \ \y/(2K ■ £) 2 - AK 2 £ 2 



(2.9) 



In another word we have 



\y/(2K -£) 2 -4K 2 £ 2 \ [ dz5{(zK 

Jo 



1. 



(2.10) 



In the same spirit of the Faddeev- Popov method in gauge-fixings, we insert it into eq (|2.3| 



d 4 £S + (£ 2 )\^/(2K ■ £) 2 - AKH 2 \ / dz5((zK - £) 2 )f(£) 



Jo 



dz / d 4 £5 + {£ 2 ){2K ■ £)5{{zK - £f)f{£). 



By use of eq ( |1.6|) and taking the QCD convention 2a ■ b = (a\b\a], then 



— c 



(2vr) 3 J 



dz / (A|dA) 



WdX 



tdtzK 2 5 [z 2 K 2 - zt(X\K\X )f(X,X,t) 

2 



dzz / (X\dX) 



(2vr) 3 J 



Putting all together we have 
d A £ 

-^(£ 2 )f(£) 



XldX 



, . , zK 2 
dtd I — - t 



K 2 



X\K\X 



X\K\X 



f(X,X,t). 



(2k) 



., , dzz / (XldX) XldX 
(2kY Jo 



, « zK2 

dtS -, r - t 



K 2 



X\K\X 



X\K\X 



f(X, X,t). (2.11) 



where c = k/2 is related to the Jacobian of changing integration variables and the way we have 
taken the residues, in which we will omit the 2ki factor. When the integrand is of the form in X4, 
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eq. ( 2. 11|) is all one needs to evaluate. When the integrand is of the form in Xe, one needs one more 



Feynman parametrization to put it in the form of Z4, as discussed in section 1.1. 

Now we make two remarks. The first is about the choice of momentum K. It is not arbitrary: 
the K must the the total energy-momentum tensor of all would-be-integrated momenta. The second 
is the integration region z 6 [0, 1]. In fact, the discussion of region of z is nothing, but the discussion 
of integration boundary faced by the evaluation of cross section. The particular simple result [0, 1] 
is special for massless particles in 4D. If it is massive, the region of z will be function of mass and 
total energy-momentum tensor K. 

2.2 Arbitrary number of out-going particles and the recursive method 

Having worked out the processes of three out-going particles, we can do the integration for arbitrary 
n out-going particles recursively 



i=l 

d 3 L n 



(2tt) 3 J 2E t 



-I n -i(f;K -L n ;L n ,Pj), (n>3). (2.12) 



where in the second line the idea is showing explicitly. More accurately, assuming / is the function 
of n should-be-integrated momenta Lj, as well as the total energy-momentum tensor K = Y^=i 
and other external momenta Pj, the integration can be done in two steps. At the first step, we 
leave the L n un-integrated, i.e., we integrate other (n — 1) momenta Lj. For this case, the total 
energy-momentum tensor should be (K — L n ). After the first step, we integrate the left L n again 
using ([2.11] ) but now with the total energy-momentum tensor K. Again, there is no need to carve 
out a complicated integration region in the phase space for massless particles. The integrations are 
always simply over the interval [0, 1]. 

The idea of recursive evaluation is very natural and has been applied in other methods, for 
example, straightforward evaluation using momentum components. This direct evaluation will be 
simplified if we choose right reference frame, for example, the center-of-mass frame. In such case, it 
usually takes efforts to rewrite results after integration as the Lorentz invariant form. This becomes 
more severe when we have multiple out-going particles. 

In contrast, when the spinor integration method is used recursively, Lorentz invariant expres- 
sions are obtained automatically at each step. There is no need to specify any frame and this 
simplifies the calculation process greatly. 

Finally we want to make some observations. First, for each J d 3 p integration, we have one 
J dz integration, thus for n out-going particles, we have (n — 2) ^-integrations (there maybe other 
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Feynman parameter integrations from spinor integrations). When (n — 2) is large, usually we can 
not find the analytic expressions, but since for each z the integration region is [0, 1], the numerical 
evaluation should be easy to realize. 

Secondly, When converting from momentum variables to spinor variables in all integrations, 
there will be Jacobi factors floating around. In addition, when taking the residues of poles, we have 
omitted the 27ri factor. That is to say, normalization factors will be needed in the conversion. In eq. 
( [2. 11| ), it is c = 7r/2. Fortunately, for given number of out-going particles, this normalization factor 
is universal. The simplest way to fix these normalization factors is to calculate the pure physical 
phase space volume by both methods, as to be shown in the following examples. 



3. Example one: two out-going particles 

Let us start with the simplest case of two out-going particles. As we have mentioned in introduction, 
it has exactly the suitable form for the spinor integration method. In the following three subsections, 
we will use the spinor method to perform the integration. A brief review of method has been 



presented in the introduction and more details can be found in references |]?fl- [[L2|| . For comparison, 
we also list the results from standard momentum integration. In the final subsection, we will briefly 
discuss the IR/collinear divergences. 

As we have mentioned in the introduction, although these examples are very simple and bear 
little practical importance, they are good for the demonstration of our method, from which the 
character of our method is clear. The validity of our method for general input shows in the Appendix 
B. 

3.1 The pure phase space integration 

We denote the physical phase space integration of n out-going particles as J* orm (/; K), where s 
stands for the spinor method and m the momentum method. The K is the sum of momenta 
of these n particles and / can be a function of K and other external variables as well as the n 
would-be-integrated momenta. 

Spinor integration method : The integration is given by 

J tdt J (\\d\) \\d\ 5((K-L) 2 ) 



—c 



(27T) 5 

where c = ir/2 is the Jacobi factor of changing integration variables. As mentioned at the end of the 
last section, we will omit the 2m factor when taking the residues of poles. It has been calculated 
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in various references J?] and the result is 



Momentum integration method: The expression is given by 

3 



/AT i 
(2^2E^ K - L >- L ^ 



(3.2) 



Taking the center-of-mass frame, where K = (E, 0, 0, 0), and by use of 
5{E -E x - E 2 ) 



2E X 

we get 



5((E - E 2 y - Ef) = 5((K -L 2 - la) ) = 5(E 2 - 2EE 2 ) 



which is identical to the result obtained from the spinor method. It is worth to notice that a suitable 
choice of reference frame has been made to simplify the calculation. 

3.2 The example with / = 2L X • L 2 

Spinor integration method: With this integrand, we have 7 

l2(u;K) = J ^y 3 S + (Ll)5 + (Ll)(2^5\K -L 2 - L X )1L X L 2 
= WyS rf4L ^ + (^)^ + ((^-^) 2 )2L 2 (K-L 2 ). 

Using L\ = there is only one term left. Following the process in the first section, we use ( |1.6| ) to 
rewrite the integration and integrate t with delta function, then use ( |1.15| ) to get the final result 



Momentum integration method: Similar to the pure phase space integration, we have 

™> = (2^ = W?2 K2 (3 ' 5) 

which is identical to /J (it; K). 

7 In fact, using 2L\ ■ = {L\ + L2) 2 = K 2 , this example is same as the pure phase space integration. Here we do 
it using different way to demonstrate the technique of spinor integration. 
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3.3 The example with / = l/(2P 1 ■ L 2 + m 2 ) 

This example is actually of certain practical value. In real cross section calculations, we meet 
infrared/collinear divergences when there are massless particles. One way to regularize such diver- 
gences is to add a mass term in the propagator. 
Spinor integration method: The integral is 

Using ( |1.6| ) to rewrite the integration and integrating t with delta function, we have 

where a 2 = m 2 /E 2 and the normalization factor 7r/2 has been inserted with the spinor integration 
variables. Introducing a Feynman parameter, we rewrite the above formula as 



^A|(l -x)K + x{P 1 + a 2 K)\\ 
The result can be written down directly 

1 7T 1 7T 1 1 1 

l2{f;K) = j2^2j dX ((1 - x)K + xJF, + a 2 K)) 2 = 2 WfK 2 " ln (1 + (3 ' 8) 

Momentum integration method: Choosing the center-of-mass frame with K = Pi + P 2 = 
Lx + L 2 , we have following components K = (E, 0, 0, 0), P x = (E/2, 0, 0, P), P 2 = (E/2,0,0,-P) 
and Li = (E\, —ki), L 2 = (E 2 , ki), thus 2P X L 2 = EE 2 (1 — cos9), 9 is the angle between Pi and L 2 . 
Doing the integration we are left with 

k ^ = ^aF sin m ^ — 1 ^vt 2 = ln (i + i)- ( 3 - 9 ) 

8{27i) J - cos 61) + m 2 (2tt) 4E Z a z 

Notice that the result has a logarithm-dependence divergence when m — > 0, which is the usual IR 
divergence related to massless particles. Again in this calculation, special reference frame has been 
chosen while in spinor method, each middle step is Lorentz invariant. 

3.4 The divergent behavior 

To discuss further the divergent behaviors of cross-sections for massless particles, we now consider 
a simple physical process, in which two gluons are scattered another two Pi + P 2 — > P3 + P4 of 
different helicity configurations. We will focus on the divergent behaviors in these calculations. 



XldX 



\\K\X 



1 

AlPi + a 2 K\X 



(3.6) 



7T 



(2vr) 2 2 



dx / (X\dX) X\dX 
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Given color structure, there are three amplitudes with different helicity configurations. Two of 
them are 

(1 3) 4 



A 12 (l-2-3 + 4 + ) 



(i 2r 



A 13 (l-2+3-4+) 



(2 3) (3 4) (4 1)' -"v- - - - , (1 2) (2 3) (3 4) (4 1)' 
First consider A 12 and denote = (Pi + Pj) 2 . In the QCD convention we have 8 



(3.10) 



C2 

2 °12 



fl2 - \A 12 \ - 2 



(3.11) 



'23 



By the momentum method, we have 

d 3 P 3 d 3 P 4 1 1 



(2tt) 3 (2tt) 3 2E 3 2E 4 



(2tt)^(K - P 3 - P 4 



(2PiP 2 ) 2 
(2P 2 P 3 ) 2 



Taking the center-of-mass frame and setting the z axis along the direction of pi, we have P\ = 
(E, 0, 0, P),P 2 — (E, 0, 0, —P), thus the integration is reduced to 



/ 2 m (/i 2 ; K) = — 



sin Ode 



(3.12) 



2tt J 2(1 - cos Of 

A singularity appears when p% is parallel to p 2 and 6 = 0. This is the familiar collinear singularity. 

The same divergent behavior can be observed in the spinor method. After some elementary 
algebra, we get the following integration form 



1 7T 



(2tt) 2 2 



(X\d\) \\d\ 



S 



12 



A|P 2 |A 



2 P 2 ' 



(3.13) 



Just as we have seen from the momentum method, a divergence appears, since P 2 = 0. 
The situation for A13 is a little more complex. We have 



13 



113 



o4 

D 13 
C2 q2 ' 
°12' ;, 23 



The momentum method gives 



2n J 32 (1 — cost') 2 



It can be rewritten as 



I?(f 13 ;K) = J- [\m9d9 ( 1 

27r J Q \2(1 — cost 



1 



+ 



1 



1 



1-cos^y 2(2tt) 2(2tt) 12(2tt) 



(3.14) 



s It is just one term of the real cross section since we have not added up all color-ordered contribution. 
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where we have left the first two infinite terms un-integrated. 
Now take the spinor method. One has 



Wis; 10 



(X\dX) 



XldX 



Sm(X\1\X 



AI2IA 



X\K\X 



/ 



7T 



(2tt) 2 2 



6S 



12 



(A|dA) X\dX 
AS l2 (X\2\X 



S 



12 



4S 



12 



A|2|A 
S 12 (X\2\X 



X\K\X 



X\K\X 



+ 



A|2|A 
2 \ 



X\K\X 



X\K\X 



J 



(3.15) 



where we have used the spinor algebra to split the first line into the second line. 

It is not difficult to see that the five terms in I 2 (fi3] K) and 7™(/ 13 ; K) correspond to each other 
exactly. In particular, the divergence of 1/(1 — cos 9) is related to the divergence of 1/ (A|2|A . 



4. Example two: three out-going particles 

4.1 The pure phase space integration 

Let us start from the pure phase space integration and identify the normalization factor c = 7r/2, 
as mentioned in section 2.1. From eqs (2.3.20) and (2.3.38) of |13|], the phase space integration is 
reduced to 

l ^ = T§^S? Xl llj x * = S«r- (4 ' 1) 

where K 2 = Q 2 . In this formula, Lorentz invariant variables x±,x 2 has been used by initial sim- 
plification. This is possible because for small number of particles, there are only a few Lorentz 
invariant quantities we can construct. With more and more particles, the number of such quantities 
will increase dramatically. 

To make use of the spinor method, we have the integration 

The / d i L 2 d i L^ + {L 2 )5 + {Lf ) (2'n)H\{K - L x ) - L 2 - L 3 ) is just the form of the phase space in- 
tegration of two outgoing particles, which we know the result is 1/87T. Putting it back we are left 
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only with 



7T 



(2vr; 



IT —C 

2 (2^ /, 
7T c/r 2 



cfez / (X\d\) X\d\ 




(4.3) 



2 2(2tt) 5 ' 

Comparing both results we find immediately c = 7r/2. 
4.2 The example with / = (Li + L3) 2 

First take the spinor method. Notice that s = (L\ + L 3 ) 2 = (K — L 2 ) 2 , which can be used to 
simplify the calculation significantly, since the integrations over L\ and L 3 are the same as that of 
the pure phase space integration. Then we have 



11 (r,K) 



d A U 



1 



■5 + (L 2 2 )(K 2 -2KL 2 ) ■ — 



57T 



-It 



2(2tt) 5 



dz 



z(l-z) 



(X\dX) 



XldX 



2\3 



7l 2 (K 



2\2 



X\K\X 



24(2tt) 



(4.4) 



Similar to equation (J4.1|), one has s = (L\ + L 3 ) 2 = Q 2 {1 — x 2 ) (as given in |ll|), thus 

Q 2 



' 16(2^ , n 

4.3 The example with / = s 2 t 

The integration is 



d%\ ( dx 2 Q 2 (l — x 2 ) — 



2\2 



l — Xl 



96(2tt) 3 ' 



(4.5) 



n{s 2 t-K)= J] ( p^S^L 2 )(2n) 4 5 4 (K-L 3 -L 2 -L 1 )(2L 1 L 3 ) 2 .2L 2 L 3 



(4.6) 



This calculation is a little more complicated. However, the basic steps are the same as those for 
/ = s. First, integrate over L\ with the delta function of the energy-momentum conservation; then, 
integrate over L 2 with the remaining delta functions in spinor coordinates; finally, integrate over 
L 3 according to eq (|2.11| ). The process is generic. 

For illustrations, here we write out the main steps 

I s 3 (s 2 t;K) = JL y d'L 3 5 + {L 2 3 ) J d*L 2 5 + (L 2 2 )5 + ((K - L 3 - L 2 ) 2 ) x 
((2KL 3 ) 2 ■ 2L 2 L 3 - 2(2KL 3 )(2L 2 L 3 ) 2 + (2L 2 L 3 ) 3 ) 
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7T 



(2tt) 5 2 

1 7T 2 



1 2 1 

2 ~ 3 + 4 



d 4 L 3 5 + (L2) ( 1-1 + 1 ) (2KL 3) 



XldX 



2\5 



X\K\X 



Using the momentum method [13|] we have 

Q 2 



I?{sH-K) 



16(2tt) , 



i ,i 

dx x I dx 2 {Q 2 ) 2 {l-x 2 ) 2 Q 2 {l-x l ) 

Jl~x 1 



7l 2 (K 



2\4 



240(2vr) 



(Q 



2\4 



960(2vr) 3 ' 



which is identical to t; K). 



(4.7) 



(4- 



5. Example three: four out-going particles 

In this section we show how to deal with multiple out-going particles recursively. To be concrete, 
we will focus on the two-four process. Generalization to cases with arbitrary n-outgoing particles 
will be straightforward. 

5.1 The pure phase space integration 

Using K) as obtained in the previous section and the recursion relation, we have 

1 7T 2 



11(1; K) = j£*sj d*L 4 5 + (Ll)I°(l;K-L 4 ) 



4(2tt) 8 2 



J d 4 L 4 5 + (Lj)(K-L 4 ) 



1 vr 2 



K 2 / d^L^+iLl) - I d 4 L 4 5 + (Ll)2KL 4 



4(2vr) 8 2 

Integrating out these two terms separately, one has 



2\2 



768(2vr) E 



(5.1) 



which can compare with the calculation using momentum integration method. Other method, like 
using the Lorentz invariant variables as = Li-Lj, can also be used to do the calculation. However, 
unlike the situation with only three L iy thing becomes more complicated. 

5.2 Two examples of non-trivial / 

The first example is f\ = s, and the second example is / 2 = (Li+L2) 2 (L2+L 3 ) 2 (L 3 +L 4 ) 2 (L 4 +Li) 2 = 
(2LiL2)(2L2L 3 )(2L 3 L 4 )(2L 4 Li), which is more complicated. Here the integrations will be performed 
with the spinor method. The momentum method can be used to make comparisons. For fx, the 
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evaluation is trivial and will not be presented. For / 2 , the calculation is a little involved and will 
be deferred to the Appendix A. 

For the fi = {L x + L 2 ) 2 , we have 

JJ(s; K ) = J^j d*L 4 5+(Ll)I°(s; K - L<) = J *L<S* {L\) , (5.2) 

by using the result for I^(s;K) obtained in the last section. The last integration is easy to do. 
Following the prescribed procedure, we get 

(K 2 ) 3 



It(s-K) 



(5.3) 



4608(2tt) 5 ' 

Now we turn to the function f 2 = (2LiL 2 )(2L 2 L 3 )(2L 3 L 4 )(2L 4 L 1 ). The factor (2tt)~ 8 and the 
Jacobi factor will be dropped for simplicity in the following. The first step is almost the same as 
that of the momentum method (notice for I 2 , the K is in fact the K — L 3 — L 4 ), 

m 2 ;K) = J d 4 L 2 6 + (L 2 2 ) J d A L 1 6 + (L 2 1 )6 4 (K-L 2 -L 1 )(2L 1 L 2 )(2L 2 L 3 )(2L 3 L 4 )(2L 4 L 1 ) 

= J d 4 L 2 5 + (Lj)5 + ((K-L 2 ) 2 ) 

x ((2L 3 L 4 )(2KL 4 )(2KL 2 )(2L 2 L 3 ) - (2L 3 L 4 )(2KL 2 )(2L 2 L 3 )(2L 2 L 4 )) . (5.4) 

The first term is simple and it is (K 2 )(2L 3 L 4 )(2KL 3 )(2KL 4 )/2. 

The second term is a little bit more complicated. To proceed, we convert the momentum 
measure into the spinor measure, and perform the t-integral with the 5-function: 

d 4 L 2 5 + (Lj)5 + ((K - L 2 ) 2 )(2L 3 L 4 )(2KL 2 )(2L 2 L 3 )(2L 2 L 4 ) 



(2L 3 L 4 ) / (\\d\) \\d\ 



{K 2 ) 4 (X\L 3 \X (X\L 4 \X 



(5.5) 



X\K\X 



Then, we do an auxiliary integration as suggested in section 1.1: 



(X\dX) X\dX 



(K 2 ) 4 (X\R\X 



X\K\X 



n 1 



(5.6) 



with R = xL 3 + yL 4 . From eq. ( |1.16| ), one has 

1 



I aux = -{(K 2 )(2KR) 2 -(K 2 ) 2 R 2 ). 



(5.7) 
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Next, pick the coefficient of the (2xy)-term in I aux , which is 



\{K 2 )(2KL 3 )(2KL 4 ) - l -(K 2 ) 2 (2L 3 L 4 ). (5.8) 

6 O 



Thus, the end result of ( p.5| ) is 



\{K 2 ){2KL 3 ){2KL 4 ){2L 3 L 4 ) - \{K 2 ) 2 {2L 3 L 4 ) 2 . (5.9) 
Combining the two terms in ( |5.4|) , we have 

I°(f;K-L 3 -L 4 ) = \{K 2 ){2KL 3 ){2KL 4 ){2L 3 L 4 ) + hK 2 ) 2 {2L 3 L 4 ) 2 . (5.10) 

o o 

There is no more difficulty in remaining steps. After some elementary manipulations, we have 
I°(f 2 ; K-U) = f d 4 L 3 6 + (Ll)I°(f 2 ; K-L 3 -L 4 ) = ^L(K 2 f(2KL 4 ) 2 , (5.11) 

im K) = j d*L 4 6 + (L 2 4 )I s 3 r (f2; K-L 4 ) = j^iK 2 ) 6 - (5-12) 

Putting all factors back, we get finally 

I'(fcK) = l -——,(K 2 Y. (5.13) 

4U ' ; 675 x 2 12 (2tt) 5V ; v ' 

as to be checked with the calculations in the Appendix A. 
6. Conclusion 

In this paper, we have proposed to use the spinor integration method, which is developed from the 
phase space integration of one- loop unitarity cut, to do the real physical phase space integration 
for the total cross section. This paper, as the initial construction of frame, focus on the massless 
particles in pure 4D. We will discuss the massive and general D-dimension case in subsequent work. 

Now, let us sum some salient points of our method. The first point is the rewriting measure 
given in ( 2.11|) by the Faddeev-Popov trick. This new form enable us to do the spinor integration 



method. Since the ^-integration is always [0, 1] for massless case, we can change the order of 
integration and leave all ^-integration at the end after we have performed all spinor integrations. 
With n out-going particles, there will be (n — 2) ^-integration to do. The integrand will be the 
Lorentz invariant expression of momenta of in-coming particles only. With two particles pi,p 2 we 
have only one nontrivial variables K 2 = (pi+P2) 2 - For massive case, the integrand for z will depend 
on K as well as various masses and thing will be more complicated. 
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In comparison, there will be (3n— 4) integrations in general, if we use the momentum integration 
method 9 . Furthermore, the expressions are usually in component-form which is not manifest 
Lorentz invariant. With more and more out-going particles, it will also become difficult to specify 
the integrated regions and separate angles variables and module variables. 

One can also try to rewrite the measure using Lorentz invariant variables. With n Li there are 
n(n 2 ~^ quantities given by = Li ■ Lj. They are not independent to each other in general, so the 
discussion of proper choice of subset as well as the constraints among them becomes more and more 
tedious with the increase of number of Lj. 

Having shown the promise of our new method, the result in this paper is not immediately to be 
useful for practical calculations. As we have mentioned several times, we need to include massive 
outgoing particles as well as general D-dimension. This can, in fact, be accomplished by mimicking 
one-loop calculations, where one generalizes pure 4D spinor integration to (4 — 2e) dimension spinor 
integration ||, [Kj] although some technical difficulties need to be attacked. Work on this is in 



progress. 
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A. Using momentum method to integrate fi for four out-going particles 

In this appendix, we will present the integration with f 2 = (2L 1 L 2 )(2L 2 L 3 )(2L 3 L i )(2L 4: Li) for four 
out-going particles by using the momentum method. It is to be compared with the spinor method 
presented in section 5.2. For simplicity, the factor (2tt)~ 8 has been dropped, just as in section 5.2. 
In the spirit of recursion, we first perform the integration 

d 3 L 2 fd 3 L 1 



2£^ J ^ 54 (^- L 2-^i)(2L 1 L 2 )(2L 2 L3)(2L 3 L 4 )(2L 4 L i ; 
= J^{{K-L 2 f) 



x ((2L 3 L 4 )(2KL 4 )(2KL 2 )(2L 2 L 3 ) - (2L 3 L 4 )(2KL 2 )(2L 2 L 3 )(2L 2 L 4 j) . (A.l) 

In the following, K will be substituted by K — L 3 — L 4 . The first term is simpler than the second. 
It is 

' ! ' L2 -6 + ((K - L 2 ) 2 )(2L 3 L 4 )(2KL A )(2KL 2 )(2L 2 L 3 ) 



2E< 



'With a suitable choice of reference framework, we may reduce dimension of integration somewhat further. 
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= (2L 3 L 4 )(2KL 4 ) J ^r5 + (E 2 - 2EE 2 ){2EE 2 )(2E 2 E 3 {1 - z)), (z = cos# 23 ) 

= (2L 3 L 4 )(2KL 4 ) • — (E 2 )(2EE 3 ) = — (K 2 )(2KL 4 )(2KL 3 )(2L 3 L 4 ). (A.2) 

o o 

In the last line, we have transformed from the component-form to the Lorentz-invariant form, such 
as 2EE 3 — > 2KL 3 , to make it suitable for the use of recursion relations. Such manipulations will 
always be used in the following. 

Next we deal with the second term (without the minus sign) . The difficulty is that the integrand 
depends on two angles, 623 and 824, which cannot be integrated over at the same time. However we 
can spilt the integral into two terms, either of which has only one angle. To do this, we introduce 
two auxiliary parameters L + = L 3 + L 4 and L_ = L 3 — L 4 . Then the second term becomes 

f d 3 l 

J ^d 6+iiK ~ L 2) 2 )(2L 3 L 4 )(2KL 2 )(L 2 (L + + L_))(L 2 (L + - L_)) 

j^S + ((K - L 2 ) 2 )(2KL 2 ) [(L 2 L + ) 2 - (L 2 L_) 2 ] . (A.3) 

Now we calculate the first term of the expression only. The second term is almost the same 
as the first. Set L + = (E + , P + p + ), L_ = (E_, P_p_), where P + and P_ are the module of the 
corresponding momentum vector. 82+ is the angle between / 2 and p+, and 82- is that between Z 2 
and p_. Then the first term becomes 

(2L 3 L 4 )^ dip j dz + j ^LdE 2 5 + (E 2 -2EE2)(2EE2)(E 2 E + -E2P + z + ) 2 , ( 2 + = cos# 2+ ) 

= (27r)(2L 3 L 4 ) {^{EE + f + ^pA . (A.4) 

Also we have L 2 + = E\ — P\ = (L 3 + L 4 ) 2 = 2L 3 L 4 , i.e. P\ = E\ — 2L 3 L 4 . Substituting it into eq 
flAll), we have 

(27T) (^(K 2 )(2L 3 L 4 )(KL + ) 2 - ±(K 2 ) 2 (2L 3 L 4 ) 2 
Similarly, the second term gives 

(2vr) (^(K 2 )(2L 3 L A )(KL_) 2 + ^(K 2 ) 2 (2L 3 L 4 ) 2 
Combining them together, we get the result of the second term in ( |A.1|) : 

(2tt) f^(K 2 )(2KL 4 )(2KL 3 )(2L 3 L 4 ) - ±(K 2 ) 2 (2L 3 L 4 ) 2 ) . (A.5) 
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Then 

I?(f2;K) = (JL-JL\ (K*)(2KL,)(2KL 3 )(2L 3 L,) + -(K 2 ) 2 (2L 3 L 4 ) 2 

= 2 ^(K")(2KL A )(2KL 3 )(2L 3 L 4 ) + |^ 2 ) 2 (2L 3 L 4 ) 2 . (A.6) 
The rest calculations are tedious but not difficult. The end results are 

-^I?(h; K - L :s ) = _i_(2,r) 2 (A' 2 ) 3 (2A'L 4 ) 2 , (A.7) 



Adding back the factor of (2tt) 8 , we finally get 

IT(f 2] K) = 1 - x A E 2 f. (A.9) 

4 kji, i 2764800(2vr) 5V ; v ' 

which is the same as /K/2; K), as obtained in section 5.2. 

B. The spinor integration as integration in complex plane 

In this section, we will rewrite spinor integration as the integration in the complex plane. The 
correctness of the spinor integration method will be obvious. Furthermore, as complex integration 
can be applied to any form of inputs in principal, so can the spinor integration be. This is very 
important since for phase space integration, the input has singularities other than those in the 
propagators in the one- loop calculation (see, for example, fl5(| ). We will demonstrate it with a few 
examples. 

B.l Rewriting of spinor integration 

Since spinor A has only two components, we can expand it with two independent spinors, for 
example, A a , as 

|A) = |a> + 2|6>, |A] = \a] + z\b) . (B.l) 

For real momentum, z is the complex conjugation of z. Using this expansion, we have (A dX) = 
(a b) dz, A dX = [a b] dz, thus the measure is given by 

J (\\d\) \\d\ = (a\b) [a 1 6] J dzdz . (B.2) 
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With an arbitrary integrand G(X, A), using above replacement we have 



(X\dX) 



XldX 



G(X, X) = (a\b) [a\b] / dzdzG(z,z) 



(B-3) 



In another word, the spinor integration is a two-dimensional integration over complex plane. 



For integration in complex plane, there is an important formula [14 



or in the form 



j) (udx + vdy) 



gdz 




dv du\ 
r \dx dyj 




dg 



dzAdz . ls - form - 1] 



ic J JR,dz 

Now we see how to apply above formula with the simplest example: 

1 /"/",,.,.,.. 1 



(B.4) 



(B.5) 



(A dX) 



A dX 




(a b) [a b] dzdz 
— (b (a b) [a b] dz 



{x\K\xy 

d 1 




(a b) [a b] dzdz 



((X\K\a] + z (X\K\b}) 2 



-1 



dz(X\K\b] ((X\K\a\ +z(X\K\b}) 
1 -1 



(X\K\b) ({X\K\a]+z(X\K\ 



(B.6) 



Using (A|iT|6] = (a|i^|6] + z t^b\K\b\ with pole location z 
-=4 (where we have neglected the 27ri factor). 



{a\K\b\ 
' (b\K\b] 



. Taking residue of pole we get 



It will be useful to compare above calculation with the one we did in spinor integration. We 



have 



(A dX) 



X dX 



J (x\K\xy 



(A dX) 



dX <9j; 



T] X 



(X\K\ V }(X\K\X 



(B.7) 



Taking the residue of (Al-fT)^] we have -j^. If we set rj = b and using 



oz 



dX\dr 



(B.8) 



we will see explicitly the one-to-one correspondence between our spinor integration method and 
complex integration method. 

The power of (p.5| ) is that it reduces the two-dimension integration into pure algebraic calcu- 
lation, i.e., reading out residue by implicitly using 

d 1 



dz z — a 



%5 2 (z — a) 



(B.9) 
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Spinor integration method has also used this property, which is casted into the spinor formula and 
called the "holomorphic anomaly" ||. 

It is obvious now that any spinor integration can be rewritten as complex plane integration, so 
its validity is obvious. However, keeping the spinor form the expression will be much more compact 
and also easier to read out residues. 



B.2 More examples 

The spinor method is developed for the one-loop calculation at beginning, thus the most cases 
we have met are propagator- like singularities, where we have result ( |1.15|) . For other kinds of 
singularities, we need similar results. This can be obtained by first going to complex variables and 
then pulling back to spinor variables. Using this technique we can do the spinor integration for 
other kinds of singularities. Let us demonstrate this with a few examples. 
First let us discuss 



(A dX) [X dX] 



(X\K\X\ (X\R\X] 



(a b) [a b] dzdz 



(a b) [a b] dzdz( 



((X\K\a] + z (X\K\b])((X\R\a] + z (X\R\b\) 
(X\K\b] (X\R\b] 



(X\K\a\ + z (X\K\b\ (X\R\a\ + z (X\R\b\ 



) 



(X\K\b\ (X\R\a] - (X\R\b] (\\K\a] 
d 



(a b) dz—(\n((X\K\a] + z (\\K\b\) - ln((A|i2|a] + z {X\R\b})) 



x 



(X\KR\X) 



(B.10) 



Using (|B.8|) as well as 



a 



d_ 

'OX 



(3\X =|)9|a] 



(B.ll) 



above result is equivalent to 



dX 



(X\K\R\X) 



In 



X\K\X 



X\R\X 



X dX 



X\K\X 



X\R\X 



(B.12) 
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Similarly we have 



/ 



dX d~ x 



{\\K\R\X) 



i 



\\R\X 



\\K\X 



X dX 



X\K\X 



X\R\X 



Using these results, we can do more examples: 

Example 1: The first example is the square root form as 



h = J d'L^+^^-Kf)^ 
= [ (X dX) 



R 



X dX 



K 2 I ~l3/2 

X\K\X 



X\R\X 



(X dX) 



V R2 



dX\-t 
. dX. 



X\R\X 



X\K\X 



Reading out residues we have 



where A = (2R ■ K) 2 - AR 2 K 2 . 



2(V2K ■ R- + V2K ■ R + y/A) 



(B.13) 



(B.14) 



(B.15) 



Example 2: 

The second example is the one with gram determinant. This is, in fact, the one we will meet in 
one-loop calculation as well as in tree-level after one external particle having been integrated out. 
It represents the typical singularity behavior in the phase space integration. 



J d 4 L5 + (L 2 )5 + ( 



)5 + ((L - Kf) = J d 4 L5 + (L 2 )5 + ((L - Kf 



( Kf Ki ■ K 2 Kx ■ K 3 \ 



det 



K x ■ K 2 K 2 K 3 ■ K 2 

-2 



\K x .K 3 K 3 .Kt K 2 J 

d 4 L5 + (L 2 )5 + ((L-K) 2 ) 

2(K -R)(L- R)(K ■ L) - R 2 (K ■ L) 2 - K 2 (L ■ R) 2 
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where we have taken K\ = K,K 2 = R, K 3 = L. Doing standard manipulation we get 

-4 1 



I 2 = J (A dX) [A dX) 



K 2 2{K ■ R) (X\K\X\ (\\R\X\ - R 2 (X\K\X} 2 - K 2 {X\R\X\ 



<" a H"*1f(wt»^ (B16) 



where 



-K-R+ JTR ■ R) 2 - K 2 R 2 

2/ 2 = l, yi = -R 2 , Xl = K-R+y/{K-R) 2 -K*R?, x 2 = V y R2 ' 



Now this is the familiar case given by (|1.13|) and we have the final result 



K 2 J ((1 - a)(arii2 + y x K) + a(x 2 R + y 2 K))' 



da— r (B.17) 



K 2 ((K ■ R) 2 — K 2 R 2 ) J a(a-l) 

which is divergent with a-integration. This example is, in fact, another typical example of divergence 
at the boundary of phase space. 

It is worth to emphasize that the factorization QB.16 ) we have done is the general method we 



need to use for phase space integration in more complicated situations. 
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